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Abstract 


Hydrostatic adjustment, due to d inbat. ic heating in two nonisot -hernial atmospheres is examined. In the 
first ease the temperature stratification is continuous; in the second case the atmosphere is composed of 
a warm, isothermal troposphere and a colder, isothermal, semi-infinitely deep stratosphere. In both cases 
hydrostatic adjustment, to a good approximation, follows the pattern found in the Lamb problem (semi- 
infinitely deep, isothermal atmosphere): Initially we have acoustic waves with the kinetic energy increasing 
or decreasing at the expense of available elastic energy. After this initial period the acoustic waves evolve 
into acoustic-gravity waves with the kinetic, available potential and available elastic energies interacting wit h 
each other. Relaxation to hydrostatic balance occurs within a few oscillations. 

Stratification in an atmosphere with a continuous temperature profile affects primarily the shape and 
amplitude of the disturbances. In the two-layer atmosphere, a certain amount of energy is trapped in the 
tropospheric waveguide as disturbances reflect off the tropopause and back into the troposphere. With each 
internal reflection a portion of this trapped energy escapes and radiates to infinity. 



1. Introduction. 


fhe atmosphere continuously undergoes various forms of adjustment. On the large scales baroclinic 
waves are excited and reduce the north-sout.fi temperature gradient so efficiently that the extratropical 
atmosphere remains near neutral stability (Stone 1078). Inertial-gravity waves radiate ageostrophic distur- 
bances out to infinity, leaving behind a geostrophic final state ( Ross by 1038). On the mesoscale hydrostatic 
imbalances are eliminated by the excitation of acoustic-gravity waves. 

Although baroclinic and geostrophic adjustment have received considerable attention, it is only recently 
that hydrostatic adjustment has been studied. Bannon (1995) examined how hydrostatic imbalances excite 
acoustic-gravity waves in an isothermal atmosphere, the so-called Lamb problem. An important aspect of 
his analysis was a computation of the energetics. Bannon showed that nonhydrostatic disturbances may 
generate three types of energy: kinetic, available potential and available elastic energies. Although he found 
the initial and final distribution of these three energies during hydrostatic adjustment, he did not follow their 
temporal evaluation. We will do this in Section 2. 

Although Lambs problem provides many useful insights into hydrostatic adjustment, the atmosphere 
is not isothermal. Vertical stratification plays an important role in determining the vertical propagation of 
acoustic-gravity waves from their origin in the troposphere into the upper atmosphere (Francis 1973. Fried- 
man 1962). Consequently, vertical stratification must also play an important role in hydrostatic adjustment. 
The purpose of this paper is to examine those effects. 

To realize our goal, I have examined two stratifications which highlight different, effects. In Section 3, 
we examine hydrostatic adjustment when the temperature profile is given by T,)(z) = T m ( 1 + This 

basic state yields a nearly constant lapse rate in the troposphere and an essentially isothermal stratosphere: a 
similar model was examined by Pekeris ( 19-18). With this thermal stratification, acoustic-gravitv waves travel 
at different phase speeds within the atmosphere, just as electric waves do within a nonuniform transmission 
line. 

To examine the case when a sharp troposphere is present, we introduce in Section 4 a t wo- 1 aver at- 
mosphere where one layer represents the troposphere and the other models the stratosphere. Both lavers 
are isothermal but at different temperatures. In this case acoustic-gravity waves will reflect at the interface 
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(tropopause) and energy leaks from the troposphere inti 1 * the stratosphere. Finally. in Sort ion 5 vve present. 


our conclusions. 


2. The Lamb Problem 


As a benchmark for later computations we first examine the energetics during hydrostatic adjustment 
in an isothermal atmosphere. The hydrostatic adjustment arises due to diabatic heating of finite vertical 
extent. 

We begin by modeling the atmosphere as a stably stratified, inviscid. diatomic, ideal gas in a Carte- 
sian coordinate system where —g k equals the acceleration due to gravity. The one-dimensional linearized 
equations for small-amplitude perturbations about the basic state 

Po(-) = pu ( - ) = p.er z,H -\ 


T 0 (r) = T.. 6 0 (z) = 


( 2 . 1 ) 


are 


dw dp 

P 0 lH = -^- piJ - 


dp dp { ) die 

P0 _ + __, + , o _ 


= 0. 


dp 
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+ P° c '~q7 ~ P°3 W = -^Pl H ( r 
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H{: — 2a)]S(t). 


( 2 . 2 ) 

(2.3) 

(2.4) 


and 


% + = 6 j ^-[H{z - a) - H{z - 2 a)]S(t). (2.5) 

at az 7Po 

where the basic state quantities have a subscript nought and the asterisk subscript denotes a constant 
reference value. Here w denotes the velocity in the vertical ^ direction at time t . p denotes pressure. I 
denotes the temperature, 6 denotes the potential temperature, p denotes the density. II 5. the scale height, 
equals RT*/g. R denotes the ideal gas constant, 7 = C P /C V = 7/5 is the ratio of the specific heats at constant 
pressure and volume, c 2 = 7 RT* and k = (7 — I ) / 7 . The functions II ( ) and S( ) denote the Heaviside step 
and delta functions, respectively. In addition to the governing equations, vve have the boundary conditions 
that 


ir(i).t) = Q and lim pt } ( z)u 2 (r. t) — > 0 


( 2 . 6 ) 
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and the initial conditions that 


p(r.O) = »(-,()) =/>(;.()) = ir(c.0) = 0. 


Equations (2.2)- (’2.5) may be combined together to yield 

0_ po w 2 ‘ P 2 po f <j0\ 

dt 2 2/>oe- 2A r ,y \0o/ 


,V 0 ==f *L = e. (2,8) 

0 00 dz Hs 

Equation (2.8) describes the energetics of the model; the first term on the left side of (2.8) gives the kinetic 
energy per unit volume while the second and third terms give the available elastic and available potential 
energies, respectively. See Bannon (1995). 

To solve (2.2)-(2.5), we combine these equations together to obtain a single equation in w: 


6 2 w dw nd~w dip , 

*^T + 91-z = d (- - «) ~<H- - 2a )] d (0- 

at- dz dz- Po 


(2.10) 


In appendix A. I solve (2.10) by Fourier sine and Laplace transforms. Once w(z.t) is known, the other 
solutions follow from (2.3)— (2.5). These solutions are 
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A7 - T 0 Ap/p., Aw - Ap/p.c, sgn( ) is the sign function and J 0 ( ) and J x ( ) denote Bessel functions of the 
first kind and zeroth and first order, respectively. The integrals in (2. 14)— (2. 15) were evaluated using Cote’s 
sixth order scheme (see Froberg 1968, p. 201). 

For a fixed value of a/H s . p(zj), w(z,t), and d(z.t) may be computed. Initially waves are excited. As 
these disturbances propagate out to infinity, the atmosphere approaches the steady-state aiven bv 
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Although we could present the temporal evolution of these fields, we have chosen to proceed along a different 
path. From (2.8). we have that 
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the next section where the temperature deceases with heigh, as T. (1 + 


t. The Exponential Temperature Profile 


Having found the waves excited m an un 

a continuously varying temperature field, namely 


iform (isothermal) atmosphere by a region of diabatic heating. 


we now examine 
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This particular temperature profile produces a constant lapse rate in the lower troposphere 

. .sotherma, aimosphere m the upper stratosphere (r » «• " 
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in at the ground to its minimum 
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a/ L) — I ■'). (•?. I ) yields ;i surface temperature of 2*<J K and temperatures of 2a9 K ami 2 i:> K ;U the 
of o km and 10 km. respectively. 

Our analysis again starts with (2.2)--(2.. r )). Again these equations may be combined together 
(2.10) with the exception that now the basic state varies according to (d.i), 
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where p. = p.RT . and p, is the pressure at c = 0. For a statically stable atmosphere. D/H s > 3/k. 
The nondimensional. Laplace- transformed vertical velocity equation becomes 
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In Appendix B 1 solve (3.5). In terms of the nondimensional r, H (r.s) is 
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if a < r < 2, and 
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if 2 < r. The parameters a', b' and c' are defined by (B.4)-(B.5). In the limit of £ -4 0 (an isothermal 
atmosphere), the hypergeometric functions tend to one and we recover Lamb's solution (A. 5). During the 
numerical inversion of (3.6)— (3.8) a numerical algorithm developed by Mavrhofer and Fischer (1994) was 
employed to compute Gauss hvpergeometric series for moderate values of |s|. For large |.s| . the asymptotic 
expansion (Luke 1969, Section 7.2) for the series was used. 

To obtain p(z,t) and 6[z,t)* we take the Laplace transform of (2.4) and (2.5). The nondimensional. 
Laplace-transformed form of these equations are 


sQ(z.s) 


(*(.- I) 2) 


and 


sP(z.s) = H{; - 1) - H{z - 2) + pr){z]c ’ ( '" s) 
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AH of the transforms were inverted numerically using Talbot's method (Murli and Rizzardi. 1990). except 
for those terms when the inverse could be obtained by inspection. 

To illustrate our solution I first found the T m for a given 3 which gave a surface temperature near 288 
K. Because /„ decreases with 3. so does // s and aj increases. I hese values are summarized in Table 1. 





lor a given .) and «/ lf s I then computed the energetics for large I. Only the available potential and 
elastic energies are nonzero and they are given in ['able 1. The case of J = 0 provid.nl a check on the 
code. As this table shows, as J increases the available potential energy decreases while the available elastic 
increases. Taken together, their sum decreases with a corresponding increase in the amount of transients. 
This is shown in Fig. 3. Although there are very slight changes in the locations of the maxima and minima, 
the evolution of the various energies follows the life cycle found in the Lamb problem. 

To understand these slight changes I have plotted the vertical velocity at various times in Fig. 4. As one 
might expect, differences in the local speed of sound has affected the propagation of the various wave fronts, 
resulting in slight changes in the exact placement of the maxima and minima in the energetics diagram. On 
the other hand, we can identity the same progression of wave pulses regardless of the value of J. 

Having found the effect that continuous stratification has on hydrostatic adjustment, let us look at the 
limiting case when the stratification is concentrated in a discontinuity - . 

4. Two-Level Model 


As counterpoint to the results presented in the previous section where the basic state's temperature 
decreased exponentially from 7.(1 + 3) to T.. we now examine hydrostatic adjustment in a two-level atmo- 
sphere consisting of an isothermal troposphere at temperature 7j (or scale height /7j) and depth H which 
lies underneath an infinitely deep, isothermal stratosphere at temperature T. (or scale height H s ). Because 
each layer is isothermal, (2.2)-(2.5) still holds in each layer. However, the basic state is now given by 
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( sing the same nondimensionalization as in Section 2. 1 show in Appendix C that the Laplace transform 
of the nondimensional velocity is 
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The variable /i z is given in ( C . / ) with the restriction that Re(/io) > 0. We may compute 
from (3.9)— (3. 10) - 


(4.7) 

© ( r , s) and P(z. s) 


Equation (4.5) (4./) provide a clear physical picture of the wave field. The first term in (4.5) represents 
the direct wave excited at z — z s . Waves radiate both upward and downward. The downward propagating 
waves eventually reflect at the earth's surface and then propagate upward. This process is given by the 
second term in (4.5). 

W hen a direct or reflected wave reaches the tropopause, a portion is reflected back into the troposphere 
[the term given by Ae~ J ^ lH / a in (4.5)] and some of the energy escapes into the stratosphere [the term given 
by .4e j n (4.6)]. The waves that are reflected by the tropopause propagate down to the earth s 

surface where they are reflected. These reflected waves again propagate to the tropopause where some of the 
energy is reflected and some escapes. Thus the terms given by A n €~‘ 2n,il Ff l a in (4.5) represents the internal 
reflections within the tropospheric waveguide while the corresponding terms in (4.6) represent the continual 
bleeding of energy out of the troposphere. Eventually all of the transients escape, leaving the hydrostatically 
adjusted basic state. Because these modes represent the leakage of wave energy away from the troposphere, 
we shall refer to them as "leaky modes'. 





To illustrate our results I have computed solutions when the tropospheric temperature is 290 k and the 
stratospheric temperature is 150 and 220 K. Although 220 K is a realistic temperature for the stratosphere, 
150 K is too cold. My reason for including t his case is to see what happens when the temperature contrast 
across the tropopause is very strong. As in the previous sections. Talbot's method was used to numerically 
invert the Laplace transform term by term. 

Figure 5 shows the temporal evolution of the kinetic (I\E). available potential (PF) and available elastic 
(AE) energies for the stratospheric temperatures of 290, 220 and 150 K. The tropospheric temperature 
is maintained at 290 K. As this figures shows, until the wave front reaches the tropopause, the temporal 
evolution is the same, regardless of stratospheric temperature. The differences in the curves of AE and PE 
are due to different values of steady-state AE and PE. Table 2 gives these steacly-state values and shows 
that they increase with a decrease in stratospheric temperature. 

The general evolution of the KE, PE and AE is very similar to the Lamb problem. However, as 
the stratospheric temperature decreases, several additional and shorter period oscillations develop: this is 
especially evident in the plot of kinetic energy. These additional undulations are due to the reflections within 
the tropospheric waveguide and the subsequent release of energy into the stratosphere on a later reflection. 

To examine how the presence of a strong tropopause affects the actual wave solutions I have graphed 
the nondimensional vertical velocities at various nondimensional times C\t/a . where c i is the speed of sound 
in the tropospheric layer. Figure 6 shows these results when the stratospheric temperature is 290 and 150 K. 
The effect of the colder stratosphere is two-fold: First, wave solutions that propagate into the stratosphere 
travel at a slower speed. This is quite reasonable because the temperature and speed of the sound are lower 
there. Second we see multiple reflections in the amplitude of the vertical velocity from reflections off the 
tropopause. 

5. Conclusions 

The effects of stratification on hydrostatic adjustment arising from tropospheric diabatic heating have 
been examined. These are the key finding: 

• To a good approximation, hydrostatic adjustment in a stratified atmosphere is similar to that in an 
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isothermal atmosphere: Initially acoustic waves are generated with the kinetic energy growing or decaying 
at the expense the available elastic energy. Within two oscillations the acoustic waves evolve into acoustic- 
gravity waves. Hydrostatic adjustment occurs during the next to I oscillations. 

• Itt an atmosphere with continuous stratification, the primary effect of the stratification is on the shape and 
amplitude of the disturbances as differences in the local speed of sound accelerate and retard propagation. 

• In an atmosphere with a sharp tropopause, small, but distinct, reflections may be observed. These internally 
reflecting waves trap a small amount of the energy within the tropospheric waveguide. Eventually all of the 
transients are able to escape and radiate out to infinity. 
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Appendix A 


We begin our solution of (2.10) by first, introducing the nondimensionnl variables r ' = z/a. /' = ct/a 
and w' — w/Au\ Dropping the primes, we have that 
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where l'F(r.s) = e a: l- Hs VV’(r, s). To solve (A. 2), we use Fourier sine transforms. We have chosen this 
particular transform because it ensures that tt’(0,/) = 0 for all f. Upon applying this transform, we obtain 
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where we have evaluated (A. 4) using integral tables (Gradshtevn and Ryzhik 1965. formula 3.742.1 ). Because 


expf-^Cs 2 + a 2 ) 
s 2 + a- 


Jo(a\/< 2 -6 2 ) H{t-b ), 6 > 0. 


(A. 6) 


we obtain (2.11) from (A. 5). 

To compute p(z.t) and 9{z.t). we take the Laplace transform of (2.4) and (2.5) and find that 
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f(») = 


0 -6v^ 5 J + el- 


and G(s) 


e -6V'-» J + a- 

$>/ 5- + a 2 


b > 0. 


To invert F{s). we note that 


f -bs _ « 2 + a- _ 


cib C 


Ji ( a 0“ — ) 

-//(/- 6) 

V 7 /" - 6 2 


or 


so that 


^ -6\/ s-'+a 2 
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//(t-6) f 

Jb 


cr 


- 6'/ ..- + <(■ 


H(t - b) 


1 -hias/v’ -b'-) 

-ab f 

Jb v/r- - V- 


b s/T^V 

' ^(aV- 2 -/>-) 


dr 


dr 


In a similar manner. 


’-1 
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/ o . 

SVS" + a- / Jb 


J Jo r 2 - 6 2 ^ c/r. 


Appendix B 


We begin our solution by introducing \\ (c,s) = e at/z ^ s \\'{z, s). Eq. (To) then becomes 


, _ a - / 0 \ / cf 2 VV 2<i ^ r/VV o “V ~ 

1 + Je ) f + — — + _VV 

ac~ //.s' fi: //“> 


a ( d\Y av ^ 

TT S \ 7T + 777 >V 


.s-VV 


Po(c 


— :/H* [«»'(- _ !) -S( : -2)}. 


We now introduce the new independent variable it = ~ . Equation (B.l) becomes 

c/->V d\V .s „ n-„ n u o / //.? - 1 

(1 " K) " + [c' - ( 1 - - 6'H— - a'-6'-VV = - «i) - ^(« - >*2) 


(A. 8) 


(A. 9) 


(A. 10) 


(A. 11) 


(A. 12) 


(A. 13) 


(B.l) 


(B 2) 



wlic'ro 


1 1 I -l.s-//; 

t/ = \/ L -f — ^ 

2 2 V a- 




■in 


+ 777'/ + 


(0.3) 

(0.4) 


(0.5) 


and Ui denote the value of it corresponding to z t . 

Our particular choice of u was motivated by two factors. First, this choice yields the classic Gauss 
hvpergeometric equation. Second, we anticipate that the solution for the exponential temperature profile 
should be very similar to Lamb’s solution. After all, the temperature profile differs markedly only near the 
ground. If we choose the v given by (B.3), then e auz i 5 becomes the Lamb solution and VV(c.s) is merely 
the correction arising from the e~ a: t D term in the temperature profile. 

Consider the differential equation 


( i \ d 2(J 

(1 - u)u — 
au- 


+ w 


(1 


b')ll 


di 

du 


2 b ,2 y = -6(u-<;). 


(B .6) 


With the requirement that #(— ,J|0 = 0. 


■ f Ci[F(a',6';c'; u)G(a ' . b ' . c': -3) - F(a\b' ,c'\ -3)G(a* ,b':c': u)], 

C\F(a'.b'-c',u). 


-3 < u < C- 

c < u < 0. 


(B.7) 


where F[a\ b':c f : u) and G(a f . b': c'\ u) are hvpergeometric functions of the first and second kind, respectively. 
See Lebedev (1972, Chapter 9). The arbitrary constants C\ and Co are found by requiring that 


lim «|C) — lhii </(w|£) and (1 — u)u-~ 
u -K“ u-k+ du 




(13.8) 


The second condition in (B.8) follows from integrating (B.6) over the infinitesimally thin interval [c . C* 
where (J - and denote points just below and above the point u — These relationships yield 


Ci - 


r(a , )r(6')F(a , 1 6':c / :0 


T(c')0- C '{l - O ci - a '- h ‘F(a'.b / :c'\-J) 


(B.9) 


and 


Co = 


r(a')r(&')[/>'. V: c':QG(a\ b’ . c': -J) - F(a'.b'. c f : -J)G(a'.b': c': £)] 

-Q c ‘- (l, - h ' F(a',b’\c f : -J) 

l l 


(B 10) 



whore L( ) is the gamma function (Lebedev, 1972. Chapter l). In deriving (B.9) (B. 10) we have used Uie 
properties of the Wrou.sk ian as they apply to hypergeomet rio functions (Lebedev, 1972. p. 278. problem 9). 
Substituting of C\ and 0> into (B.7) yields 


(-J) l " c ' F(a\ 6': o': u)F{ l + a' - c\ 1 + V - c\ 2 - c'; -J)F(a'.b': e'X) 

~ (c' - 1 )C 1_L “' ( 1 - C) 

^-^(f/.frV'i-cQF^ + a! - F. L -f F - F;2 - F: u)F(F. F: F: Q 
(F - l)C l - c ' (l — F(u'. F; F: —J) 

if — J < u < £, and 

(-j) l - c/ F(a',^;c':()F(l + a' - F, 1 + F - F,2 - F; -3)F(a\b'\ F: u) 
~ (F - l)C 1 ~ c/ (l ~;)c'-a'-*'F(a',F;F:-J) 

C 1 ~ c F(a / .F > c / : - J)F(1 + a' - F, 1 + V - F; 2 - F; QF(a'. b': F: a) 

(F - I)C 1_C '(1 -C ) c '- a '- b ' F{a',b':c':-3) 


(B 11) 


(B 11) 


if C < « < 0. To derive (B. 1 1)-(B. 12) we used the relationships that r(F)/F(F — 1) = c' — 1 and the 
relationship between G(a\ F, F; c) and F{a\ 6\ F; r) (Lebedev. 1972, p. 277. problem 7). I he final solution 
is obtained by substituting for a', b f and F, applying the results to each of the impulsive forcings in (B/2) 
and multiplying by the appropriate coefficients. 

Our ability to express G(ci\ F:c r) in terms of F(a'. F; F; z) and F(1 + a f — F. 1 + F — F; 2 — F; r) 

requires that the contour of integration not. pass near the singularities F — 0. ±1, ±2 When the contour 

generated by the numerical inversion did, we took an equivalent contour [increased A as suggested by Murli 
and Rizzardi (1990)] to avoid coding this special case. 


Appendix C 


We begin our analysis by taking the Laplace transform of (2.10). Because of the linearity of the problem, 
we only have to solve an equation of the form 


H s d 2 w x 

H x OA 


a dw i 

Jh~d7 



u " / S(z - z s )S(t) for 0 <z<H/a 
Hipn(z s ) 


(C.l) 


and 


d 2 W 2 t a dw< 2 cFlt'j _ n 

i)A + 7Ts~0z 0A~ = 


for 


H/a < z. 


b'> 


(C.2) 



To find t. in; solution for (’2.10) we tin'll subtract, tin? solution for r, = 2 front the one for - 1. Again I 
have scaled the vertical distance with a and time with <t/c„. Denoting the transform of the nondimensional 
vertical velocity by I4*(c..s). we have that 


^ VV, = r----*/-"' ~M= ~ ->,) 

d:- \ /f, 4 Hipo('s) 


for the troposphere and 


, 3 + __ VV . = ° 


for the stratosphere, where = e a ^~ Hl VV'i (c, s) and lb>(-^') — e a ^- W “ s H 2 (-^). 

Taking wi(OJ) = 0 and the radiation boundary condition into account, the solutions to (C.3)-(C.4) are 


A si nh ( ! r ) , 


0<z<z s 


Vi, / \ _ I .ion ~ — — * 

1 \ B sinh(^i z) + Ccosh(/fi z). z s < z < FI/ 


VV 3 (~. s) — Dexp[— ^*j(r - FI /a )], 


•> Us *> , *> *> - 

« = 77T S + ifff a " 5 = '■ + 5«r 


Here we have chosen the branch Re(^ 3 ) > 0. The remaining task is to find .4, C and D . The kinematic 
boundary condition leads to VV*i ( //, s) = VV 3 (//. s), where VV 3 (c. s) = c #/-//* VV 3 (~, $)• A combination 

of the dynamic and kinematic boundary conditions applied to the Laplace transformed (2.4) yields the second 


condition that 


dWx(H'S) dWi{H.s) fa a 
Tz “ dz + 1 2 H s ~ 


. s). 


At the point of singularity, continuity of the solution requires that VV i ( - 3 . s) — bijr/.s). where 
and .-+ denote points just below and above c = r s . respectively. Finally, an integration of (C.3) over the 


interval [z s , z+] yields the final condition that 


_ p -az t /2H 


, FIsp. 

HiP^s) 


Substituting (C.o)-(C.6) into these four conditions gives the solutions listed in Section 4. except in terms 
of hyperbolic cosines and sines. Finally we express the hyperbolic functions in terms of exponentials and 
expand the denominator as a geometric series in 


If) 
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Tables 


TabU' l: Distribution of available potential (PE) and available elastic (AK) energies in the final equilibrium 
state in an nonisothermal atmosphere where the temperature behaves according to (3.1) with D/ II. s = lo. 
[ he percentages give the ratio of the corresponding energy to the total energy. 


3 

T. 

T 0 ( 0) 

a/ FIs 

PE 

AE 

0.0 

290 K 

290 K 

0.500 

64.74% 

6.7 1% 

0.1 

260 K 

288 K 

0.558 

60.58%' 

8.32% 

0.2 

240 K 

286 K 

0.604 

56.92% 

9.79%: 

0.3 

220 K 

286 K 

0.659 

53.57% 

11,39%. 


Table 2: Same as Table 1 except for 

a two-level atmosphere with 

a/H s = 0.5. 


Tropospheric 

Stratospheric 

PE 

AE 

Temperature 

Temperature 



290 K 

290 K 

64.74% 

6.71% 

290 K 

220 K 

65.55%' 

7.35% 

290 K 

150 K 

66.69% 

8.39% 
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Tables 


'Fable L: Distribution of available potential (PE) and available elastic ( A E ) energies in the final equilibrium 
state in an nonisothermal atmosphere where the temperature behaves according to (3.1) with D/Hs ~ 1.5. 
The percentages give the ratio of the corresponding energy to the total energy. 


3 

T. 

To( 0) 

n///s 

PE 

AE 

0.0 

290 K 

290 K 

0.300 

64.71% 

6.71% 

0.1 

260 K 

288 K 

0.558 

60.58% 

8.32% 

0.2 

240 K 

286 I\ 

0.604 

56.92% 

9.79% 

0.3 

220 K 

286 K 

0.659 

53.57% 

11.39% 


Table 2: Same as Table l except 

for a two-level atmosphere 

with a / H s — 0.5. 


Tropospheric 

Stratospheric 

PE 

AE 

Temperature 

Temperature 



290 K 

290 K 

64.74% 

6.71% 

290 K 

220 K 

65.55% 

7.35% 

290 K 

150 K 

66.69% 

8.39% 
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Figures 


Figure 1 : The transient portion of available elast ic (AE). available potential (PL) and kinetic (KL) energies 
in terms of the percentage of total energy as a function of nondimensional time ct/a. Lhe solid 
line gives the case when a/Hs — , 0.25 while the dashed line corresponds to aj H 5 = 0 . . 

Figure 2: The variation of integrand of (2.10) corresponding to the the available elastic' (AL). available 
potential (PE) and kinetic (KE) energies as a function of nondimensional distance z/a at the 
nondimensional times of ct/a = 10 (solid line), 25 (dashed line) and 50 (dot-dashed line). The 
parameter a/Hs — 0.5. 

Figure 3: Same as Figure 1 except for an atmosphere with a basic state temperature given by (3.1). The 
solid line gives the isothermal case when T1 = 290 K and 3 — 0 while the dashed and long dashed 
lines correspond to XI = 250 I\. 3 — 0.15 and T m = 220 K. 3 = 0.3. respectively. 

Figure 4: The nondimensional vertical velocity in an atmosphere where the basic state temperature is given 
by (3.1). The solid line gives the isothermal case of X. = 290 K and 3 = 0 while the dashed line 
corresponds to X, — 220 K and 3 = 0.3. The corresponding a/Hs are 0.5 and 0.66, respectively. 
The snapshots are at the nondimensional times ct/a of (a) 1.0, (b) 1.5. (c) 2.0 and (d) 2.5. where 
c is the speed of sound corresponding to 290 K. 

Figure 5: Same as Figure 1 except for a two-level atmosphere, lhe solid line corresponds to the case when 
the stratospheric temperature is 290 K: the dashed line, 220 K: and the dash-dotted line. 150 K 
In ail cases the tropospheric temperature is 290 K. The nondimensional time is now c\tja and 
a/H$ = 0 . 5 . where c 1 is the speed of sound in the troposphere. 

Figure 6: Same as Figure 4 except for a two-level atmosphere where the tropospheric temperature is 290 K 
and the stratospheric temperature is either 290 K (solid line) or 150 K (dashed line). The snapshots 
are at the nondimensional times c\tja of (a) 1.0. (b) 1.5, (r) 2.0 and (d) 2.5 and a/Hs = 0.5, 
where c\ is the speed of sound in the troposphere. The tropopause lies at z/a — 2.827. 
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Figure 1: The transient portion of available elastic (AE). available potential (PE) and kinetic (KL) energies 
in terms of the percentage of total energy as a function of nondimensional time ct/a. The solid line gives the 
case when a/H$ — 0.25 while the clashed line corresponds to a / H $ — 0.5. 


20 






NONDIMENSIONAL DISTANCE 

Figure 2: The variation of integrand of (2.19) corresponding to the the available elastic (AE), available 
potential (PE) and kinetic (KE) energies as a function of nondimensional distance zja at the nondimensional 
times of ct / a = 10 (solid line). 25 (dashed line) and 50 (dot-dashed line). The parameter a/ fls = 0.5. 
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NONDIMENSIONAL TIME 


Figure 3: Same as Figure 1 except for an atmosphere with a basic state temperature given by (3.1). I he 
solid line gives the isothermal case when T m = 290 I\ and 3 — 0 while the dashed and long dashed lines 
correspond to 7. = 250 K. 3 = 0.15 and T . = 220 K. 3 = 0.3. respectively. 
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0 1 2 3 4 5 

NONDIMENSIONAL DISTANCE 


Figure 4: The nondimensional vertical velocity in an atmosphere where the basic state temperature is given 
by (3.1). The solid line gives the isothermal case of T. = 290 K and 3 = 0 while the dashed line corresponds 
to 7T = 220 K and 3 = 0.3. The corresponding a/H$ are 0.5 and 0.00. respectively. The snapshots are at the 
nondimensional times ct / a of (a) 1.0. (b) 1.5. (c) 2.0 and (d) 2.5. where c is the speed of sound corresponding 


to 290 K 








Figure 5: Same as Figure 1 except for a two-level atmosphere. I he solid line corresponds to the case when 
the stratospheric temperature is 290 K; the dashed line, 220 K: and the dash-dotted line. 150 K. In all cases 
the tropospheric temperature is 290 K. The nondimensional time is now c\t/a and c if Hs = 0.5. where Ci is 
the speed of sound in the troposphere. 
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Figure 6: Same as Figure 4 except for a two-level atmosphere where the tropospheric temperature is 290 K 
and the stratospheric temperature is either 290 K (solid line) or 150 K (dashed line). The snapshots are at 
the nondimensional times c\tja of (a) 1.0. (b) 1.5. (c) 2.0 and (d) 2.5 and nj — 0.5. where C\ is the speed 
of sound in the troposphere. The tropopause lies at z/a = 2.827. 
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Popular Summary 


The atmosphere is full of waves which differ widely in length and time between peaks. A special grow]), 
called acoustic-gravity waves, are generated whenever large pressure differences occur between adjacent 
chunks of air, such as during earthquakes, developing thunderstorms, and even during the launch of the 
space shuttle. Although a theory has been constructed of how these waves would propagate out to space 
if the atmosphere had the same temperature everywhere, this is not very realistic because temperature 
decreases with altitude. In this paper realistic* temperatures are taken into account and it is shown how these 
variations affect this radiation of energy into space. With the advent of nonhydrostatic models which include 
these acoustic-gravity waves, this study will hopefully provide tools to better understand these new models. 


